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Abstract. This work illustrates how several new pricing formulas for exotic 
options can be derived within a Levy framework by employing a unique pricing 
expression. Many existing pricing formulas of the traditional Gaussian model 
are obtained as a by-product. 



1. Introduction 

In the last decades the hterature on option pricing under Levy processes mush- 
roomed. The reason is that Levy processes have the flexibihty to capture empirical 
characteristics in stock returns such as jumps in price due to market shocks and 
distributional properties such as skewness and semi-heavy tails, while maintaining 
most of the analytical tractability of the Gaussian model. In this paper we adopt 
the class of regular Levy processes of exponential type (RLFE) as the driving pro- 
cesses, following [7]. As [7] points out, it is the most tractable subclass of Levy 
process from the analytical point of view if the Brownian motion is not available. 
Their characteristic exponents ■0(C) ^njoy very favourable properties as symbols of 
pseudo differential operators, since their real part behaves as c 1^1*^ as \£,\ — > oo in 
the strip of regularity, with positive c and v. Thus the integrals appearing in the 
pricing formulas are absolutely convergent thanks to the terms of the form e~'^'^^^\ 
Moreover, one can differentiate under the integral sign or shift the line of integration 
by using the Cauchy theorem for holomorphic functions. Such properties allow for 
a great flexibility of the method when working out the analytical pricing formulas 
for several exotic options. Moreover, such a class incorporates most of the models 
which have been proposed as an alternative to the traditional geometric Brownian 
motion, including the Normal Inverse Gaussian (NIG) (Barndorff- Nielsen [6]), the 
Hyperbolic (H) (Eberlein and Keller [13]) and the more general Generalized Hy- 
perbolic (GH) (Eberlein and Prause [15]), the four-parameter distribution named 
CGMY after the names of Carr, Geman, Madan and Yor [10] and generalized in 
[11], just to quote the most popular ones. The development of pricing models re- 
placing the traditional underlying source of randomness, the Brownian motion, by 
a Levy process has fostered a good deal of work on exotic option pricing which par- 
allels the existing results in the Gaussian framework. (See [27], for a compendium 
of recent research on the topic). 

In this article the focus is on the pricing of European exotics and the aim is to 
present a valuation formula which is the most comprehensive as possible, in that 
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several types of options can be priced directly and no specific method has to be 
devised for each of them. The formula is tailored to valuate discretely monitored 
options, which are the most popular ones in view of the regulatory issues and the 
trading practice. However, the continuous counterpart can be derived in some cases 
(see Example 2). The method employed in this work is to start from (multi-period) 
digital options as building-blocks and then to price the more complex options in 
terms of such elementary contracts. The idea that a broad class of financial deriva- 
tives can be evaluated in terms of elementary contracts such as digital options has 
been applied in the traditional Gaussian framework (see [9] and [23]) but it has not 
yet fully exploited in a non-Gaussian modeling. By providing the non-Gaussian 
counterpart of this view, we are able to obtain new pricing formulas in the Levy 
environment and to throw a new insight into some known pricing expressions. This 
method has been anticipated in [3] . This work employs a more general setting and 
provides further examples of exotic option prices. Since a Brownian motion is a 
RLPE of order 2 and any exponential type - and thus is captured in this framework 
- each example includes the known pricing expressions of the classical Gaussian 
modeling. An emphasis is put on deriving the existing formulas of the traditional 
Black-Scholes modeling from our more general framework, so that the reader is 
confronted with the more familiar expressions. To our surprise, the Lemma which 
is employed to the purpose is not foimd in the Levy literature and therefore a proof 
is given. Some applications of the main result are given in Section 4. They serve as 
an illustration of the flexibility of our valuation formula and integrate the several 
examples developed by the author in a previous work [3] . 

This paper is organized as follows. Section 2 outlines the main definitions and 
notation concerning Levy processes and provides a Lemma linking the Gaussian 
multivariate distribution to the Fourier transform-based Levy setting. This Lemma 
will be used extensively throughout the paper and has a per se mathematical in- 
terest. Section 3 presents the general valuation formula, while some examples of 
exotic option prices are provided in Section 4. The choice of the examples has been 
dictated by the wish to provide new valuation formulas and presenting them in a 
general form (for example, N-fold compound options instead of just 2-fold ones, 
flexible Asian options, both call and put options in one expression, etc.). However 
they are by no means exhaustive with respect to the potentialities of the approach. 

2. Preliminaries 

Let us consider a Levy market, i.e. a model of financial market with a determin- 
istic saving account e'"*, r > 0, and a stock following a stochastic process St = e^^ , 
where {Xt)t>o is a Levy process. As usual in a Levy setting, the Levy process 
replaces the Brownian motion which is employed in the classical modeling of stock 
prices. Here we assume that the stock price St is e'^*, where {Xt)t>o is a one- 
dimensional RLPE, i.e. a regular Levy process of order v S]0,2] and exponential 
type [A_, A+] , A_ < < A+. This means that {Xt)t>o has a characteristic expo- 
nent V(0 which admits a representation of the form: 

(2.1) = + 

where (f) is holomorphic in the strip Im^ e]A_, A+f, continuous up to the boundary 
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ofthestrip, and(/)(0 = C|er + 0(|^r')for|^|^oo and <C{l + \^{"^) 

in Im^ G [A_,A+] with vi,V2 < v. (See [B-L] for a comprehensive and Finance- 
oriented theory of RLPE processes). In order to price contingent claims on the 
stock, it is convenient to consider an equivalent martingale measure (EMM) Q which 
makes the discounted price process e~^*'St a martingale. Let ^/^p (respectively iJjq) 
denote the characteristic exponent with respect to the historic measure P (an EMM 
Q, respectively), i.e. Ep{e^^^') = e^*'^^'^?) {EQ{e'^^') = e"*'^«(«\ respectively). 
We assume that the discounted stock price e~'^*St is a martingale under Q, that 
is, we assume that the equivalent martingale measure condition (EMM-condition) 
r + TpQ^—i) = holds. Furthermore, the additional condition A_ < —1, which is 
usually assumed to let —i belong to the strip of regularity of tp, is supposed to hold 
true. In the sequel further restrictions on X± will be posed, depending on the kind 
of exotic option under study. We recall that Q is not unique in a general Levy 
setting. We do not dive into the problem of choosing a martingale measure and 
refer to the well-established literature on Levy processes in Finance dealing with this 
topic extensively (sec [17], for example). A popular method is based on Esscher 
transform and the expression for V'qIC) is obtained by first solving the equation 
-ipp^—ih) — ipp{—ih — i) = r for h, and then letting ipqi^) — 'ipp{^ — ih) — -tpp^—ih). 
In what follows we assume that an EMM Q is chosen so that Xt is a RLPE under 
it and we will omit the subscript Q both in E and ip. 

Let g{XT) denote the terminal payoff of an option on Si at the expiry date T. 
Then the no-arbitrage price of the option at the current time t {t <T) is given by: 

(2.2) F{Su t) = E[e-^^^-'^g{XT) \ X^ = ln(5t)] 

An explicit expression for F{St. t) can be obtained in terms of the Fourier transform 
g{£,) in the complex plane. Indeed, if e'^^g{x) G L^(Ili) for some oj €]A_, A+[ , then 
the Fourier transform 'g{^) of g can be defined on Im^ = w and one obtains: 

(2.3) FiSt,t) = ^ e'^^nS,-iT-t){r+^{0)^^^)^ 

In the following sections this set-up is extended to the more general case of path- 
dependent options where the payoff g depends on a set of fixed asset price monitor- 
ing times, Ti < ... < Tm < T, so that several types of exotic options can be priced 
in this framework. 

Finally we recall that a Brownian motion is a RLPE of order 2 and any expo- 
nential type. Therefore the pricing formulas of the traditional Gaussian approach 
can be obtained in our framework for several exotics. To the purpose we need the 
following: 

Lemma 1. Let C = (ckj) denote an N x N correlation matrix. Then for any 
ojk > 0, for any real number dk and with Wk = ±1, k = 1,...N, the following 
identity holds: 

= [ n Wk]NN{widu-WNdN;WCW) 

k=l,...,N 

where denotes the N-variate multinormal cumulative distribution function and 
W is the N X N diagonal matrix with entries Wk, k = 1, ...N. 
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Proof. Consider the case Wk = 1, k = l,...N, at first. Let = (wi,...,wjv) and 
let xn(a;i,...,a;jv) denote exp(- X] u!jXj)I[_di,+oo){xi)—Ii-dN,+oc,){xN) which 

j=l,...,N 

belongs to ^^(M^). Its Fourier transform in (^f , ..,^^) e is: 

exp[? X dj^j]/ n ) with = ^j^—iwj. On the other hand, exp(—5S^C 

j=l,...,N j = l,...,N 

is the Fourier transform of (p{X) = (27r)M/Vdetc 6xp(— ^X'^C"^^). Thus the left- 
hand side term can be written as a convolution of xq and 4>, that is, as 

(2.)»/^Vd5tc /«" M-Y)e-''^'' cxp(-ir^c-ir)dy, 

which is /^^ ... /_^^^ H^gp2_g=Z)dy^...rfy^ ^ iVjv(di, ...dw; C). The general case 
follows changing to variables = "^j^j ^^"-^ integrals of the left-hand term and 
replacing each dk with dk = Wkhk and C with WCW in the previous argument. 



3. Set-up and main result 

In this section an almost universal option pricing formula is proved within the 
Levy framework described in the previous section. In particular we derive the 
arbitrage-free price for generalized multi-period exotic power digital options. Such 
options are the building blocks for a broad class of exotic options with a single 
underlying asset, because several exotic options can be expressed as static portfolios 
of these multi-period power digitals. The idea of pricing several exotic options 
by a single universal formula is due to [9], where the classical Gaussian case was 
studied. Thus this Section provides the generalization to the Levy framework and 
the resulting formula is a most comprehensive one, both in terms of the kind of 
option and of the stochastic process driving the underlying asset. Since the aim is 
to establish a uniiying set-up, the formal notation is a bit involved and is laid down 
along the lines of [9] . 

Assume that the payoff of an option depends on M fixed asset price monitoring 
times, Ti < ... < Tm, where, for simplicity's sake, Tm = T, the expiry date of the 
option. Let T denote the set of times [t,Ti, ...,Tm], where t is the current time, 
t < Ti. Let Sk denote the price of the underlying asset at the monitoring time 
Tfc and let S denote the M-dimensional vector assembling all the components Sk 
which are relevant for the option under study. If 7 is M-dimensional , then S"*" 
denotes Sj^ ...S]f and 7 is referred to as the payoff index vector. We assume that 
St = e^*, where is a Levy process. Thus, in most cases, it will be convenient to 
work directly with Xk, the value of X at T^, and with the vector X= {Xi, Xm)- 
Since we want to treat call and put options together, we introduce W, a diagonal 
matrix with all the diagonal entries Wu equal to ±1. To simplify notation, wa will 
be denoted by Wi. Then both the indicator functions I[Ki,+oa) {Xi) and I(-oc,Ki]iXi) 
may be encompassed in a unique notation, li{wiXi > WiKi), depending on the sign 
of Wi. Let K denote the exercise price vector (of dimension N) and let lAf(Y > K) 

N 

denote the iV— dimensional indicator function W > Ki). In order to give 

a greater flexibility to the approach, the exercise condition matrix is introduced. 
Such term will denote any Nx M matrix A = (a„fc) involved in the payoff, where A'' 
is the exercise dimension. For example, the payoff of discrete mean Asian options 
depends on '^Ilk=i,...,M Sk ^ K or, equivalently, Yl j^^Xfe ^ In iiT = iiT. Thus 

fe=l,...M 
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their payoff can be expressed in terms of li(wiAX> wiK), wliere A = [jj, -p-] 
is a 1 X M matrix. On the other hand, the payoff of a discretely monitored digital 
option 1m{WX > WK) can be expressed in the general form Ij^iWAK > WK) by 
taking A as the M x M identity matrix. All the useful parameters are summarized 
in the payoff parameter set P= [(7^, ■■■Jm)7 

Let i^(S't, i;T,P) denote the value of a multi-period binary whose payoff is spec- 
ified in terms of T and P and where St denotes the value of the underlying 
asset at the current time t. More specifically, the expiry T payoff function is 

M 

exp[ E 7fe^fc]-ljv(T^^X > VFK). In most applications 7^ = for A; ^ M. The fol- 

fc=i 

lowing Proposition gives the arbitrage-free current value of this generalized multi- 
period option under the assumption that the value St of the underlying asset St 
is e"^*, where Xt is a RLPE of order v €]0, 2] and exponential type [A_, A+]. The 
main result of the paper is the following: 

Proposition 1. The following valuation formula holds: 

M 

F{St,t;T,P) = -I^^^St''"' n Wn f+— ... f+~-™i-i 1 

n=l 11 

n = l,...N 

N M 
n=l k=l 

M N M M 

where ...^^) = E {Tj-Tj_i)^{ E E ankU-i E Ik) withTo = i, w„ > 

j=l n=l k=j k=j 



M 

forn=l,...N and E 

k=j 



N 

E WnOJnank + Ik 
.n=l 



e]-A+,-A_[ forj^l,...,M. 



Before proving Proposition 1 we give a pricing formula for the simple case of a 
power digital option, which is slightly more general than Proposition 1 in [3] . 

Lemma 2. Let F{St,t) denote the current arbitrage- free price of the power option 
with expiry T payoff function S^li{waXT > wK), a ^ Q, w = ±1, Xt = InSt- 
Then for any w > such that awLO + 76] — A+, — A_[ , one can write 

(3.1) F{St,t) = ^ gi|[alnS,-K]-(T-t)[r+V>(a€-i7)]irf^ 

27ri J— 00— iwu) C 

Proof. The Fourier transform of the payoff function g{XT) = e''-^ li{waXT > 
wK) is g{r}) = exp[(7 — irj)^] with Imrj = —wlj^ for any such that 

7 — wcoj ^ whenever aw ^0 and wco^ e] — A+, — A_[. Then: 

F(S, t) = ■^■^snja) r+co-iww^ ^HMSt)- ^:]-{T-t){r+il,{r,)) 

which gives (3.1) changing to variables r] + i'y = a£, and letting uj^ — wy = aco. 

Remark 1. The expression (3.1) can be rewritten in terms of pseudo differential 
operators as follows: 

F{St,t) = ,«57 cxp[-(T -t){r + i^(aD^ - i7))]/(-)(ln f ) 
where l'""' denotes the indicator function I^'^\x) = -f[o,+oo)('*^2;). Here the notation 
P{Dx), with Dx = —idx , denotes a pseudo differential operator whose symbol is 
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P(^). Alternatively, F{St,t) = f{Xt,t) can be viewed as a solution to the pseudo 
differential equation: [dt — (r + ^p{Dx))]f{Xt,t) = with the final condition 
f{XT,T) = e-f^-^liiwaXr > wK). 



Lemma 3. Let F{St,t) denote the current arbitrage- free price of the power option 
with expiry T payoff function S^lN{wnanXT > WnKn, n = 1, N), a„ ^ 0, Wn = 

N 

±1, Xt = In St- Then, for any a;„ > such that Yl Wnt^nO-nk + 7 S] — A+, — A_[, 

n=l 

one can write 



-r(T-t) N n+OO — iWNl^N r+oo — iwiuii 

(3-2) nSut) = ^^^S^ UwJ ... 



00 — lW\UJ\ 

n=l,...N 



n c 



N 

exp[i E ^niO'nXt - Kn) - ^'(t,^!, ■ • ■5^)]^^ 1 ■ -C^C AT 

71=1 

N 

where ...^;v) = {T - t)^{ E - 

Proof. The result follows by writing the Fourier transform of the payoff as a 
convolution of A'' terms and by arguing as in Lemma 1. 

Proof of Proposition 1. Let us first prove the case N = 1. 

Let P= [{■j^,...-/j^j),K,w,{ai,...aM)] with w = ±1. For m, = 1,...,M let = 
K-J:T=i<^A, l*m = T.k=ilkXk {Kl = K,j*^ = 0). Let fm{Xt,t) solve 

dtfm - {r + ^{Dx))fm = for i e [r„_i,T^], with fm{Tm,Xm) = fm+l{Tm,Xm) 

for TO < M, and fmiTM^XM) = e'^'^^'^+'^*'^li{waMXM > wKlj). In view of 
Lemma 2 and Remark 1 one has: 

f^iXi. t) = ^e"'"^'+^" r+°°-™'-'e»«[''MXt-ifM]-(TM-t)[r+V(oM£-i7M)]i(7f 
Jivi\ t) J 2-Ki J — oo—iwu ^ ^* 

Then one can prove recursively that for any to: 

M 



fm{Xt,t) = -ir^ J+-_.™%xp[z^( E a,A,-i^;;)-vi/„(<,0]|rf^ 

k=Tn 

M MM 

where ^m{t,0 = E (Tj - Tj_i)V(^ E «fe - « E 7fc) with T„_i = i, w > and 

j=m k=j k=j 

M 

E [wLuak + 7fe] €] — A-|_,— A_[ for j = m,...,M. Thus h = 1 yields the result 

k=] 

in the case N = 1. Finally, the general case is proved along the same lines, by 
employing Lemma 3 and arguing recursively. 



In the sequel the pricing formulas for some exotics arc obtained as an application 
of Proposition 1. Note that our pricing formula are new in the Levy environment. 
Moreover, our approach casts a new light also on the Gaussian case that is obtained 
as a by-product. 
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4. Examples 

In this section we provide some examples to illustrate the flexibility of the main 
formula. Most of the examples we give are new in the Levy framework. Note that 
we have chosen not to include in this section the textbook examples and the most 
elementary examples concerning popular options: some are given in [3]. A few 
lines in each example are devoted to show how the well-known pricing expressions 
of the traditional Gaussian setting can be obtained as special cases of ours quite 
straightforwardly. Finally note that the references we give here are not exhaustive, 
because the focus of this paper is on the analytical formulas, while the amount of 
work on the purely computational issues is generally omitted. 

1) Forward- start options. Forward-start options are options that are only effective 
at some pre-specified time after they have been bought or sold. They are the 
building-blocks of more complex options. For example, a cliquet option is a series 
of consecutive forward-start options, where each option becomes active when the 
former expires. A valuation of forward-start options is given in [29] in the classical 
Gaussian setting. Herewith the generalization to a Levy setting is presented. 

Let T be the expiration time and let Ti < T be the time at which the option 
starts. The strike price is set to be the underlying asset price at Ti. Thus the final 
payoff is max {w{St2 ~ ■S'ti ), 0} with w ~±.\. 

As an application of Proposition 1 we give the current value of the forward-start 
option F{Si, t) in the above-mentioned Levy setting, under the assumptions A-(- > 
and A_ < —1. 

t) = ^^s^i!^-::::: i cmht^ - T,m - m 

where wi e]0, -A_[ , 002 e]0, -A_ - 1[ if w = 1 and oji e]0, A+[ , W2 e]0, A+ -|- 1[ if 
w = —1. 

Note that the greeks A (which is equal to F{St,t)/ St) and F = do not depend 
on Si. In the Gaussian case, -0(0 = — + , the formula above becomes: 

^^*(/r-™7: 1 -p[z(. + 4)(r2 - T,)e - ^2^en 

which is wSt[Nlw{^ + f)V72-Ti) - e-^^^^-^'^N{w{^ - f)VT2-Ti), i.e. the 
formula found in [29] . 

2) Asian options. Asian options under Levy processes have been priced in a 
number of papers (see [5], [18], [21]). In this subsection we show how a pricing 
formula for geometric Asian options is easily obtained from our general result. At 
first we consider discrete Asian options - whose payoff depends on a discrete average 
of the asset price at monitoring times, Ti < ... < Tm ■ The continuous average 
case is obtained as a limit. Consider a forward-start fixed strike Asian option with 
strike price K. Note that a pricing expression for floating Asian options is easily 
obtained in view of the symmetry relationship proved in [14]. Let T = Tm be the 
maturity date and let T' = Ti be the time at which the averaging starts. The payoff 
is max(w[SM — 0), where Em = ( 11 St^)"^ and w is the binary indicator. 

j=l,...,M 

In terms of Xi — In(S't) the payoff is: 

w n e^^^^li(w;AX> w;ln(ii'))-u;i^li(w;y4X> i(;ln(Js:)),with^= [^,...,^]. 

j=l,...,M 
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Thus Proposition 1 applies with = 1 and yields the following valuation formula, 
after some algebraic manipulation: 

FiS„t) = - ^-"I""" ^ In § - E {Tj - Tj-.m^^M 

with To = t, u) A_[ (]0, A+[) if w = 1 {w = —1). Note that the analogous 

expression obtained in [18], (10), is derived throughout a different argument, that 
is considering the distribution of In(SM)- 

The pricing formula for the continuous-time monitoring case, where the geomet- 
ric average is exp[ jj^^, J j,, \n{St)dt], follows from the discrete pricing formula just 
letting M — >• 00. Note that the limit can be computed under the integral sign in 
view of the nice behavior of tp- In particular, for the continuous case, one has: 

Fisut) = -^^^ /r:r: i^iw ^-pk ^-^-ii ^(^(^ - y))dyn- 

Let us now see that our formula collapses to the know valuation formula for 
discretely monitored Asian options in the Gaussian case (see [29]). Let h denote 
the averaging frequency, that is, Tj = T — (M — j)h, j = 1, M. Then, in the 

M 

Gaussian case, E (Tj - T,-OV(€^) = -iKr - + h^e^-^^^W^^ 

M M 

because E (^ - j) = ^^^^ and E (M - = ^'(^^-'^^^-^ . Thus 
nSut) = '^'^I^~eM^kl^§H{r-^){T-r))-iT-T')^^^^ 

for any a; > 0. Splitting the integral into two integrals and changing variables, one 
gets: 

F(S, t) = ^''^'^ f+'""™" exDfOTD+ - '""tj2ii^ ge--(^-') f +°°-™" exDfzr?Li- - 

where = [In f + i(r - ^){T ~ T')]/ [a^T-^ ^f^), 

D+ = D-+ aVT^^igi and /? = r(r - 1) + (r + a^^ - M^))^, which 
yields the known formula (see [29]) by application of Lemma 1. 

Finally we point out that Proposition 1 straightforwardly applies to the more 
general flexible geometric Asian options, where the flexible geometric average is 
n ^ith 9j = ^(j)/ E^li ^(i) ^-nd ^ any non- negative function (see [29] 

for the Gaussian case). The following expression is obtained: 

nSt,t) = - ^-"X""" W+i) exp[iCln # - E (Tj - Tj_,m E 0,M 

with To = w e]l,-A_[ (]0,A+[) if «; = 1 (u; = -1). 

3) Discrete lookback options. Lookback options are path-dependent options whose 
payoff depends on the extremal price of the underlying asset over the life of the 
option. We focus on the realistic case of finite sampling lookback, that is, the 
asset price is monitored at particular dates. Moreover, we confine ourselves to 
fixed strikes options, that is, the payoff is the maximum difference between the 
optimal price and the strike price K which is determined at inception. According to 
[12] "perhaps, in exponential Levy model closed-form formulas are not, in general. 
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available for pricing these options". In the traditional Gaussian setting discrete 
lookback options have been priced in [20] exactly, and in [8] by an adjustment of 
the continuous case. 

Let Ti < T2 < ...Tm be the monitoring times and assume that the lookback 
period [Ti.Tjv] starts in the future, i.e. the current time t < Ti, and ends at 
expiration, i.e. Tm = T, the expiration date. However our method applies also to 
backward starting options. The payoff is max(5'Ti , Stm ^I^) — K for the call and 
K — min(5'T^ , Stm ) K) for the put, that is, max(wS'r^ , wStm > wK) — wK for 
any option, where w = ±1. Following [20] one can write the payoff in our notation 
as follows: 

we^'^llM-iiwA^P^li. > 0)-1m(w^^^'X > wB^p^)] - wK{1-1m{wX < wK)) 

p=i 
where 

1 if j = P 

A^P^ is an (M — 1) x M matrix whose entries are = ^ —1 if i = j ^ P 

O otherwise 

if j=P,-i¥=j 

is an M X M matrix whose entries are A^l' = —1 if i = j ■ 

otherwise 

B^P^ is an M— dimension vector whose entries are with the exception of the p*^ 
that is —InK; 

K is an M— dimension vector whose entries are all — In Jf. 

Thus Proposition 1 yields the following expression for the current value F{St,t) of 
the lookback option: 

(4.1) n,^S,e-^^^M-^^ g .+oo-i^c.i _ _ .+oo-n„.,, i j_ ^ ^ _ ^ 

V / (27rz)^^ ^ J — 00—iwuJi J—oo — iwujM 11 f ^ ^ J ^/ 

p=l 11 7=1 

rx=l,...JV 

M 

^( E - 0] - E {T, - T,_0^(- E inmi-di^-idi^+i-d^M 

n<j j=P+l n>j 

w"S«e~''*^-'"~*' v-v f+oo-iwui r+oo—iwuM 1 r_ \^ (rp _ q-< \ 

(27ri)'>^-i 2^ J-oo-iwwi ■■■ J-oo-iwuM TT ^ 2^ 3 -^j-l) 

p=l 11 i = \ 

n=l,...JV 

M 

i^iE^n-i)- E (Tj - T,-i)V(- E inMi-d^p-id^p+i-d^M 

n<j i=P+l n>j 



- ^"^"fr'w""' E f+~-.'""^ ... f+<— ™ 1 exp[-iC,lnf 

^~ n = l,...N 

P M 

-E{T,-T,_,Mj:^n-ip-^)- E m-i^-iM- EOlrf^i-d^M 



(27ri)^ J— oo— Jiowi J— oo— ituwM C 

n = l,...Af ^~ 

p M 

Em -T,-_i)V(- E UHi-d^M -wKe-^^'^--') 
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M 

with some positive w„ such that such that Un < min(A+, — A_ — 1). 

Let us see how our general formula collapses to (21) of [20]. For any p = 1, M 
denote: 

9± = [Inf + (r ± ^){ T - t)]/{ a.jT-^t), ff± = (J ± 1)^/7^^ for n = 
1, 1, hpn = (§ — ^)-\/T„ — Tp for n = M. Split any integral in the first 

sum into two integrals, the former in dS,-^...d^p_-^ and the latter in d^^j^-^-.-d^j^^, and 
then change to variables ^^cr-^/l^^^T^ = ^„ {^'n^^'^n — Tp = respectively). 
Then the integrals in the first sum of (4.1) becomes the product of: 

p— 1 



rnKn 

and 

M 

r+oo-JU)Wp+i r+oo-tU't^M 1 r ir^ / v _ 1 f 2 _ 

(27ri)*'-P J-cx>-i«;Wp+i ■- J-oo-iwuM TT f ^ \'''in"'pn 2^n 

n=p+l,...M ^ 

Employing Lemma 1 each term in the first sum is transformed into: 

where the correlation matrix Ap has typical element J Tp-xZ , tn < n, and the 

correlation matrix 0p has typical element J j," _rf , m > n. A similar treatment 
on the terms in the second sum turns each term into: 

wSte-^^'^'^-'^^'>Np{wg^-^, ...,wg^p_-^,-wq+; A*p)NM-p{whp^p+i, ...,whpM;&p) 
where the correlation matrix A* has typical element p'^^ = J t'-t" ioT n > m, 
m = 1, ...,p — 1, p* = \ ^ _ " for n < p, and p*„ = 1 for any n. Finally, the 



pn ^ V "^p-* 

integrals in the third sum are turned into 'wKe~^^'^'^~*^NMi—wqi, —wq'j^; A) 
where the correlation matrix has typical element pui. = ""^'^''7 . 

4) Chooser options 

A chooser option gives its holder the right to decide at a prespecified time (choice 
date = Ti) before the maturity T whether he/she would like the option to be a call 
or a put option. As a straightforward application of Proposition 1 (with A = /) we 
give a valuation formula for simple chooser options, i.e. the call and the put have 
the same strike price K and maturity date T. Note that the decision whether the 
option is a call or a put depends on the value of: 

Max{C{ST^ ■,K,T), P{St, ; K, T)) = C{St, ■,K,T) + MaxiKe-"-^^-^'^ - St, , 0). 
In other words the choice is: 

Call ^ Ke-^<^-^^'> < St, ; Put ^ Ke-^<^-^^^ > St,. 

The payoff can be expressed as: 
{St-K)12{St > K,St, > Ke-''^'^-'^^)) + (K-ST,)U(ST < K,St, < i^e-''(^-^i)) 
Then, in view of Proposition 1, with = 1 and M = 2, the price for the simple 
chooser option at time t <Ti can be written in the form: 

F{St, t) = A^ - K * A2 + K * As - Ai 
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where the following choice are made in Proposition 1 for each term: 



K= {Ke-'-^^-^'\K) 7 =(0,1) «;2) = (1, 1) 

K={Ker-(^'^^\K) 7 =(0,0) (wuw^) = (1,1) 

K={Ke--(^-^^\K) 7 =(0,0) {wuW2) = (-1,-1) 

K= {Ke-(^-^-\K) 7 = (0, 1) {wuW2) = (-1,-1). 



for Ai 
for A2 
for A3 
for A4 
Then 

^ (2Tri)-^ J—oo—iuJ2 J — oo — iuii ^1^2 ^ 

In view of the residue theorem Ai becomes: 

^ ZTTi J— 00 — za;2 ^2 



27ri 

On the other hand 



which under the residue theorern is transformed into: 

I e-'-t'^-" r+oo-iu-i ,4 i„(^+^(T-Ti)-(Ti-t)^(4i)J_j^ 
^ 27rj J-00-iwi «i ?1 



Thus the formula is simplified because the double integrals cancel out and the final 

expression is: 

+oo+i<^i e*fi M§+'-(T-Ti)-(Ti-t)V'(«i-i)J_WfJ_ 

J-OC + Jti,'i {1 ^IJ 

Ke--'-^-*^ r r+°°+»"2 In ^-(T-t)V(?2) J_rf^^ + 

27ri U— oo+iwg €2 

f e'«i ln(#+r(T-Ti)-(Ti-t)V(£i) 1 1 
J— 00— iwi 4i 

In the Gaussian case the formula becomes 

F{St,t) = S[N{d+) - N{-d+^j,J] - Ke-'-^'^-*^N{d:^) - Ni-d^^^J] 
with 4 = (HSt/K) + {r± ^){T - t))/{aVT^) 
dT.T, = 4i + r{T - Ti)/{aVTV^) 

which is the price for a chooser option obtained by Rubinstein in 1991 (See also 
[29]). 

5) Multicompound options. A notable example of application of our Proposition 
is the valuation of multicompound options (of order N), that is options whose 
underlying asset is a multicompound option of order TV— 1. A closed form expression 
for this kind of options was developed by [19] in the usual Gaussian setting. An 
extension to a more general Levy environment has been provided in [3] for the 
compound call options of order 2. Here we give a valuation expression which holds 
for any order. Note that multicompound options have been introduced to face the 
problem of pricing defaultable coupon bonds ([19], [4]) ; however they have a wide 
range of applications to all opportunities having a sequential nature (see [1], [2] for 
applications in real option analysis). 

Let t <Ti <T2 < ...Tjv and let FNiSt,t;Ti,Ki,wi; ...■,Tn, Kn,wn) denote the 
current value of a European TV- fold compound option expiring at time Ti, with 
strike price Ki and with as underlying asset a European (N — l)-fold compound 
option expiring at time T2, with strike price K2, ....until the final underlying asset, 
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a European option on a stock, with exercise date and price given by Tn and Kn. 
As usually Wj = ±1 represents the call/put attribute of each option. Assume 
that the stochastic process followed by the underlying stock is a RLPE. Suppose 
that is the solution to Ftq^j{S* ^Tj]Tj+i,Kjj^i,Wj+\; ...;Tm,Km,wn) = Kj for 
j = 1,...,N—1. We will comment on existence and uniqueness of Sj later on. Then 
the payoff of the multicompound option can be written as: 

N N j 

n WnSTf^lNiWNSTN > WnS%, ...,WiSti > WiS^) - J2 tl WnKjlj{WjSTj > 
n=l j=l n=l 

WjS*,...,WiSTi > WiS*) 

where S'^ = therefore our method applies. The current value of each term is 
obtained by straightforward application of Proposition 1 and the resulting expres- 
sion reads: 

{2771)'^ J— OO — J —OO — IWNOJN II ^ 

n=l,...N "~ 

. exp[-*jv(t, ^1, ••, - i)]d^i--dCN- 

^ {2mp J-oo-iwiuii ■■■ J-oo-iw,Lj, n ^ ?n"^S*J 

j=l 1 1 /=n n=l 

.exp[-*,(t,ei,..,e,)]rfei-'-4, 

3 

for some positive ojn such that J2 '^n^n s] — A+, — A_ — 1[ Vj = 1, A'', and where 

n=l 

*i(*,Ci,-,c,)= tiTk-Tk-imtu- 

k=l n=l 

Note that, differentiating under the integral sign, one can prove that dhjF^ = 
for hj = J^, j = 1,...,N. Thus dstF^ is just the first integral of Fn divided by 
St- Therefore uniqueness of S* is guaranteed for any j. Existence holds for the 
multicompound call options, while in the general case it holds only for suitable 

strike prices. 

In the Gaussian case, if we change variables ^^cr-^Tj — t = rjj , let ujj be any 
positive number and denote (In ^ + (r ± ^){Tj — t))/{ay/Tj — t) by d^, and 



^' * by p if j < k, then our formula collapses into: 



N No 

St n WnNN{wNd'^,...,wid^;EN) - II WnKjNj{'WjdJ , ...,widf;Ej) 

n=l j=l 71=1 

where is a j x j— matrix with typical elements pji^WkWj when j < k. 

Remark. (Put-call parity for compound options). Note that, starting from our 
formula, one can also verify the put-call parity relationship for compound options 
as a nice exercise. Indeed 
F2{St,t;l,K,,T,;w2,K2,T2) = 

_ Ste"*''^^-*) r+co-iw2U>2 r+cx>-wi In ll+ig^ In |l^-^(t,fi,f2-i) 1 ^£ _ 

(27ri)^ J —oo—iw2UJ2 J —oo—iuji ^1^2 
Jf2e-''<^2-t) „+(x)-ra2W2 r+oo-»wi i£iln|i-+i|2ln-||--*(t.ii>?2) _ 

(27I-i)2 J-<X-iW2U>2 J-OO-iWi ^ {^{^ S^^^ 

27r J — oc—iu)i i^i ^1* 

Let US shift the line of integration Im^j^ = oji up. Since we cross the pole at = 0, 
the residue theorem gives: 
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F2{St,t;l,Ki,Ti;w2,K2,T2) = 

(2ttT)^ J—oo — iw2UJ2 J— oo+«wi 



(2;;^p J-oo-iw2UJ2 J-oc+iuJi ^ Jjj2"Sl"?2" 



J — CO — iW2UJ2 J — oo+iwi 

27r J-oc+iuii ^ i4i"?l^ 

' I- (ZTTj)'' J — CXD — ltU2W2 ?2 

j^2e-''<'^^-*' r+oo-i«)2W2 j«2ln -R|-1'(*,0,e2) J_^f _ EieZl}IlJlp-(Ti-l)ii(a)l _ 
J —ca — iw2UJ2 €2 2iri ' J 

F2iSt,t;-l,Ki,Ti;w2,K2,T2) + FiiSt,t;w2,K2,T2) - Kie-^^^'-'\ 

In summary, wc have verified that the vahie of a eaif on an option equals the value 
of the corresponding put on the same option plus the value of the underlying option 
diminished by Kie~^'^'^^~*\ where Ki and Ti are the strike price and the maturity 
date of the compound option, which is the put-call parity for compound options. 

6) Discrete barrier options. Most analytical pricing formulas for barrier options 
assume continuous monitoring of the barrier, which corrisponds to some practical 
cases (e.g. in FX markets). However in practice the barrier might be monitored only 
at discrete points in time (e.g., at the close of the market). A discrete barrier option 
is either knocked in or knocked out if the price of the underlying asset is across the 
barrier at the time it is monitored. In the Gaussian case the pricing formulas have 
been studied by [8], [9], [16]; in the Levy process models an interesting survey is 
presented in [22], where the novel method of [16] is also discussed. (See also [24] 
for a numerical approach). In this subsection we derive a valuation formula for a 
discrete barrier option as a further straightforward application of Proposition 2. 
While there exists eight barrier options types, depending on the barrier knocking 
in or out, on the barrier being above or below the initial value of spot (up or down) 
and on the call/put attribute, we confine ourselves to a down-and-out call without 
rebate. The other cases can be treated similarly. 

Let B denote the level of the barrier and suppose that the underlying asset is 
monitored at times Tj, j = 1, M — 1 before the option expiry Tm- The payoff is 
{Stm - K)\M{STi > B,j = 1, M-1; Stm > K). Therefore Proposition 2 with 
N = M and A = 7, the M x M identity matrix, yields: 

ptq j.\ _ Jt'e"''^"^""'^ r+oo-iWN r+oo-iui 



fe=l,...M-l 

M-1 M M 

. exp[i E In § + In # - E (Tj - T,-i)V( E ^kMi-d^M 

i=i i=i k=j 

M 

with To = t, COM €]1, — A_[ , LOj > and E '^j < — A_. 



5. Conclusion 

This work introduces a comprehensive option pricing formula for a very general 

family of payoffs, which includes many market-relevant option payoffs as special 
cases. The proof is based on Fourier methods and on the theory of pseudo differen- 
tial operators that have been successfully applied in the literature for option pricing 
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in Levy models. However, one does not need to be equipped with such mathemati- 
cal sophistication in order to apply the main formula and thus the result may be of 
interest also for practitioners. The unifying formula we provide encompasses many 
existing option pricing expressions and is a powerful tool for generating new valu- 
ation expressions without effort. Wc have chosen to focus on discretely monitored 
options, as these have received little attention in the literature, despite their pop- 
ularity in the trading practice. Each example is complemented with its Gaussian 
counterpart and thus, while introducing new formulas in the general Levy setting, 
the paper may also serve as a review on discretely monitored options in the tra- 
ditional Black-Scholes setting. Finally we stress that the analytical method based 
on pseudo differential operators and integration in the complex plane generates a 
new numerical method (integration-along-cut method) which often performs better 
than the Fast Fourier Transform (see [7]). Therefore numerical computation will 
take advantage of our explicit solutions from many a point of view. 
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